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Abstract — The Na + and Ag + forms of ZSM-12 zeolite microporous were used as effective adsorbents for hydrocarbons trap. 
Diffusion of a single adsorbed component including ethylene as a light molecule or toluene a heavier one inside a narrow 
one -dimensional channel was studied. Experimental desorption profiles for these two gases desorbing from both Ag -ZSM-12 
and Na-ZSM-12 were obtained at very high heating rates comparable with the one observed during the cold start of a 
combustion engine. Numerical simulation based on Fick’s law and a transient diffusion boundary value yielded a good 
fitting of these experimental results. For both ZSM-12 forms considered in this study, diffusion activation energies E 0J and E e 
for toluene are much higher than those for ethylene. The values of these parameters are higher for the Ag + exchanged zeolite 
than for Na-ZSM-12 and related to the presence of strong Fewis acid sites in the former material. In addition each activation 
energy for Na-ZSM-12 and Ag-ZSM- 12 decreases at increasing heating rate. This was explained by thermal vibrations of the 
zeolite lattice which allow higher dijfusivity at high temperature. 
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I. Introduction 

Faced with increasingly strict emission standards, automobile and catalyst manufacturers have investigated a variety of 
approaches to reduce cold-start emissions. Attention has been focused on the start-up phase, when 70-80% of the 
hydrocarbons generated (HCs) by the combustion engine are emitted [1]. Cold-start is referring to the short time period of 1- 
2 min, after engine ignition and before the three-way catalyst reaches its light-off temperature which is around 200°C. 
Zeolites have been found to be the preferred molecular sieves and the best candidates for this application due to their stability 
under a variety of conditions and their affinity to adsorb HCs. In literature, different types of zeolites have been proposed and 
tested as adsorbents for HC emissions [2]— [1 1]. It has been reported that ZSM-12, which possesses one dimensional channels 
with 12-member ring pores, performs better in the envisaged application of hydrocarbon trap compared to more typical three- 
dimensional zeolites [9]— [14] . With the growing use of one dimensional zeolites in separation and reaction processes, the 
need to better understand the diffusion and transport properties of molecules within the channels of these materials is 
important not only from a fundamental, but also from a practical point of view. The mechanism used to describe the 
simultaneous passage of ethylene (light molecules) and toluene (heavy molecules) inside narrow one dimensional channels 
was designated as the single-file diffusion[15]. This phenomenon has become a topic of current research interest. Heavy 
aromatic molecules are strongly adsorbed in the zeolite, therefore they need a higher temperature in order to desorb from 
channels. On the other hand, light hydrocarbon molecules are more weakly adsorbed and easily desorb from channels at 
lower temperature. Thus in simultaneous diffusion, the desorption of lighter molecules such as ethylene from one 
dimensional narrow channels may be retarded by co-adsorbed heavier ones such as toluene. 

Despite the large number of theoretical studies published on the single-file diffusion over the years, experimental and 
numerical single component loading (small molecules or aromatic molecules) are still rare. It is important to understand the 
transport mechanism of these molecules within the pores of zeolites in order to design more efficient applications for these 
materials. The aim of this paper is to provide a comparison between experimental and simulation results for each single kind 
of gas after developing a model for diffusion in zeolite pores. 

II. Experimental details 

The detailed synthesis procedure and physicochemical characterization results for Na-ZSM-12 and Ag-ZSM- 12 samples with 
molar ratio Si/Al = 52 used in this research were described in our previous study [15]. Fig. 1 shows the scanning electron 
micrographs of Na-ZSM-12 and Ag-ZSM- 12 zeolites as synthetized samples. They have uniform shape and particle size; 
crystals have spherical form with average crystal diameter of 1pm related to the use of methyltriethylammonium chloride as 
organic template. This is agreement with the X-ray diffraction patterns of these hydrocarbon adsorbents that were recorded 
using a Siemens D5000 diffractometer using the CuKa (k = 1.54 184 A) as shown in Fig. 2. 
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FlG.l. Scanning electron micrograph of synthesized sample (A) Na-ZSM-12 and (B) Ag-ZSM-12 (52). 



Fig.2. X-ray powder diffraction patterns of as-synthesized Ag-ZSM-12 and Na-ZSM-12 samples. 


The high intensity of the XRD lines and the low background intensity point out that the samples are highly crystalline. The 
BET surface area and micropore volume were determined using an Omnisorp-lOO automatic analyzer after degassing the 
sample at 150 °C for at least 4h under vacuum (10' 5 to 10' 4 Torr). Using Fourier transform infrared spectroscopy of 
chemisorbed pyridine, Na-ZSM-12 has a higher measured C Br 6nsted/CLewis ratio compared to Ag-ZSM-12 indicating larger 
content of Lewis acid sites. Bronsted and Lewis acid centers are detected respectively by peaks around 1550 cm' 1 and 
1448cm' 1 (Fig. 3, Table 1). The acid site density of zeolites was quantified by ammonia thermodesorption. All experimental 
results are reproduced here from our previous work (Table 2) [15]. 
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Table 1 

FTIR RESULTS OF CHEMISORBED PYRIDINE 


Sample 

Si/Al 

Peak area (a. u.) 1448 cm' 1 

Peak area (a. u.) 1550 cm' 1 

Ch rousted / CL ew j s 

Na-ZSM-12 

52 

0.012 

0.322 

35.5 

Ag-ZSM-12 

52 

0.482 

0.152 

0.42 


Cgrtinsted/CLewis — 1.88*A 15 5 o/ 1.42*Ai 448 23 


Table 2 
NH 3 TPD results. 


Sample 

Tpeakl (°C) 

NH 3 (mol/g cat 
*10 5 ) 

Tpeak2 

(°C) 

NH 3 (mol/g cat 
*10 5 ) 

TpeakS (°C) 

NH 3 (mol/g cat 
*10 5 ) 

Na-ZSM-12 (52) 

92 

32 

150 

13.4 

345 

204.6 

Ag-ZSM-12 (52) 

145 

38.5 

213 

131.4 

360 

64.7 



Wavenumber (cm-1) 



Wavenumber (cm 1 ) 

FIG.3 FTIR SPECTRA OF CHEMISORBED PYRIDINE (DESORPTION TEMPERATURE 150°C) OF SAMPLES WITH 

Si/Al=52; (A) NA-ZSM-12; (B) Ag-ZSM-12. 


Here we report the application of laser diode thermal desorption (S-960 LDTD) sample introduction source, coupled with an 
atmospheric pressure chemical ionization (APCI) source mounted on a quantum ultra AM triple quadripole mass 
spectrometer for analyte detection. This device allows reaching ramping rates similar to those reached in a muffler which is 
around 3-5°C/s. Samples preliminarily loaded with adsorbed single -component ethylene and toluene separately at room 
temperature were deposited in the Laze well plates of the LDTD source (Phytronix). 

III. Single-Component Gas Desorption Modeling 
3.1 Model description 

The desorption process of a single -component gas can be simulated using a transient diffusion boundary value problem 
coupled with Fick’s law [16], [17]. In the present study, due to the structure of Na-ZSM-12and Ag-ZSM-12samples, a one- 
dimensional diffusion problem is considered [9], [15], [18]. For this purpose, let Land t T denote the zeolite sample length 
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and the total time of the desorption process, respectively. At a given time t < 


oj. 


’ T, 


, let /?(J,x) denotes the gas density 


at the position X e [0, L] . The transient diffusion problem is given by the following equations [17], [19]: 

— ’ ^ +J/v (■/(?, x)) = 0 , V(f,x)e 0,f r Jx[0,L] 


dt 

J (t,x) = -®(t,x)V/?(f,x) 

with the Neumann and Dirichlet boundary conditions: 

J (x = 0, ? ) = 0 , V? 

[/?(x = L,f) = 0, Vf 


and the initial condition: 


/?(x,f = 0) = /? 0 (x), Vx e [0, L] 


( 1 ) 

(2) 

(3) 

(4) 


In expressions (l)-(2), (li\'(.) and V are divergence and gradient operators, respectively. j[t,X ) and Oj(t, x) represent 

the gas flux and the diffusion coefficient. Substituting the flux expression (2) into the formula (1) leads to the following 
partial differential equation: 


^ ^ = div(co(i,x)y p(t,x)} , V(f,x)e 0,f r Jx[0,L] 


(5) 


In the sequel, the material is assumed to be homogeneous. Therefore, the diffusion coefficient is only time-dependent and it 
is defined according to the Arrhenius law [16], [17], [20]: 


co(t) = CT exp 

v kT 


(6) 


where ^(fc/.mor 1 ), k[j .mol' 1 .K~ l ) and T(K) are the activation energy, the Boltzmann constant and the 
temperature. T and C are two constants. F is defined as follows [17], [18]: 


T = exp 


f E A 


\ kT f J 


(7) 


Moreover, during the experiments the temperature is increased according to a constant heating rate : 

r(V) = Tq + A t 


( 8 ) 


where T {) is the initial temperature. 

In some research works [16], [17], [19], [21], it has been shown, through numerical simulations, that the above mentioned 
diffusion model fails to predict the general shape of the gas flux, especially at the beginning of the desorption process. In 
order to overcome such a drawback, it has been suggested in the work of Chatterjee el al. [17]— [19]to use two diffusion 
coefficients: i) the first one defined inside the physical domain and ii) an additional diffusion coefficient defined on the 
domain boundary with a higher activation energy. Thus, two diffusion coefficients will be considered: 
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G)(t) = C a T a exp 
H(?) = C E r H exp 
r w = exp 


f 


\ 


V kTj 
f E^ 

V kTy 


, Vxe[0,L[ 
, if x = L 


f E„ A 


r E = exp 
E^ > E„ 


\ kT f J 
f E- A 


(9) 


\ kT f ) 


For more details about this extra-diffusion barrier the reader is referred to the following references [16], [17], [20]. 

3.2 Time and Space Discretization 

The boundary value problem (l)-(5) is solved using the finite difference method [22]. To this end, the spatial domain [0, L] 


0 ,t T 


is also discretized into 


is discretized into a grid of nodes (x. ) . such that X 1 = 0 and X NS = L . The time domain 

/ \NT 

a set of time points \tj). Q • Without loss of generality, nodes and time points are supposed equidistant with fixed space step 
Ax and fixed time step At . Therefore, one can show that: 


\tj=jxAt, Vj = NT 

[x ( . = Xj + (/ - 1) x Ax, Vz = 1, NS 


(10) 


Let denotes by pf the approximation of the solution p [ t - , X ( j at given time t . and point A' ( . Using a central difference 


scheme, the right hand term of equation (5) is approximated by the following scheme [22]: 




( 11 ) 


For the time discretization of equation (5) an explicit scheme is used [22] Thus, the time partial derivative is approximated as 
follows: 




pr'-pi 


dt At 

Combining expressions (11)-(12) leads to the following finite difference approximation scheme: 

At 


p t j+1 = p{ +0){t j )—^(p 1 M -2 p/ + pi , ) 

(Axj 


( 12 ) 


(13) 


Therefore, knowing the solution at time step t. , the solution is updated in a straightforward way for the time step t j+l . In 

expression (13) the time step At is considered small enough compared to the square of the space step Ax in order to fulfill 
the stability condition for the explicit scheme [22] . 
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IV. Experimental and Numerical Results 

4.1 Experimental results 

In our previous study using Fourier transform infrared spectroscopy of chemisorbed pyridine characterization, it has been 
shown that Ag-ZSM-12 has low measured C Br 6nsted/CLewis indicating largest content of Lewis acid sites. On the other hand, 
Na-ZSM-12 has a lower 1448 cm' 1 peak area, a characteristic band of Lewis acid centers, than Ag-ZSM-12. Consequently 
the ration C Br6 n S ted/CLewis fo r Na-ZSM-12 increased. The Ag exchanged zeolite has the highest content in Lewis acid sites and 
the lowest content in Bronsted acid sites (lowest 1550 cm' 1 peak area)[15]. 

Three high heating rates have been used for the desorption process: 3°C/s, 5°C/s and 9°C/s. Fig.4 reports TPD profiles of the 
single-component, either ethylene and toluene, from Na-ZSM-12 and Ag-ZSM-12. All desorption profiles of both ethylene 
and toluene display a single peak, except for the desorption curves of ethylene from Ag-ZSM-12 (Fig. 5), which show two 
peaks. In the latter case the low temperature peak is assigned to the ethylene desorption from Bronsted acid sites whereas the 
high temperature desorption peak is from Lewis acid centers. 




T (°C) 


Fig.4. Single component (ethylene or toluene) TPD profiles of (A) Na-ZSM-12 and (B) Ag-ZSM- 

12 SAMPLES AT DIFFERENT HEATING RATES. 
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Fig.5. Comparison of experimental and numerical data of single-component TPD profiles in 
AG-ZSM-12: (A) ETHYLENE AND (B) TOLUENE AND NA-ZSM-12: (A’) ETHYLENE AND (B’> TOLUENE WITH 

F F 

3°C/S AS HEATING RATE. VALUES OF w AND H REPORTED IN TABLE 4 WERE USED WITH 

= ex P (^J and r ? = exp (i %) 

According to Table 3 the amounts of both ethylene and toluene detected are essentially not affected by the heating rate which 
is an indication of results reproducibility. However a clear dependence of desorption temperature on heating rate for toluene 
and ethylene was observed for the two zeolites. At increasing ramp rate from 3 to 9°C, the maximum desorption peak 
temperature is raised by 10-25°C. Regarding the desorption temperature of both ethylene and toluene molecules, it seems that 
the desorption temperature of Na-ZSM-12 is always lower than the desorption temperature of Ag-ZSM-12. This result 
emphasizes the importance of the strong Lewis acid sites in Ag-ZSM-12 and confirms that the zeolite exchanged with Ag + is 
a better candidate as hydrocarbon adsorbent as indicated by its higher desorption temperatures [9], [15]. 


Table 3 

Single component TPD results 


Gas 

Rate (°C/s) 

Desorbed amounts (mmol/g) 

Desorption temperature (°C) 

Na-ZSM-12 

Ag-ZSM-12 

Na-ZSM-12 

Ag-ZSM-12 

Ethylene 

3 

0.021 

0.042 

63 

(70*; 105**) 

5 

0.023 

0.044 

73 

(81*;114**) 

9 

0.020 

0.042 

80 

(89’;126”) 

Toluene 

3 

1.73 

1.94 

112 

126 

5 

1.71 

1.92 

126 

137 

9 

1.72 

1.95 

133 

144 


Desorption temperature obtained from the low T peak 
Desorption temperature obtained from the high T peak 
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Independently of the zeolite and heating rate, the amount of toluene adsorbed and the desorption temperature are much 
higher than those of ethylene. This indicates that part of ethylene is either not adsorbed at all or so weakly adsorbed that it is 
desorbed at room temperature. It is found that toluene is more strongly adsorbed than ethylene on both Bronsted and Lewis 
acid sites [15]. As shown in Table3, the amount of ethylene adsorbed on Ag-ZSM-12 is higher than the amount adsorbed on 
Na-ZSM-12 which also indicates that the silver form is a better adsorbent for the smaller gas molecules. 

4.2 Numerical results 

In order to show the ability of the transient diffusion boundary value problem (l)-(2) to predict the obtained experimental 
results, numerical simulations were carried out using the finite difference method (13). As mentioned earlier, two diffusion 
coefficients will be considered for each gas (9). The ZSM-12 zeolite has a typical channel length nearly equal to 5 |um and a 

o o o 

pore size of 5. 7 x 6.1 A . Molecular diameters of ethylene and toluene are 1.4 A and 5.7 A , respectively. 

In Fig. 5 (A’ and B’), numerical simulation results for ethylene and toluene desorption process using Na-ZSM-12 zeolite are 
compared to experimental trends. It is shown that the diffusion model predicts desorption curves in a good agreement with 
experimental results. Values of the activation energy used for these simulations are reported in Table 4. For the three heating 
rates, activation energies associated to the ethylene are smaller than those of toluene. This trend could be explained through 
the dimension of ethylene molecule which is smaller compared to the dimension of zeolite channels, in contrast to the 
toluene molecule. As a consequence, ethylene can desorb with a lower energy of activation. For both gases and for all 

heating rates the activation energy inside the domain (E^) is lower than the activation energy at the lattice 

boundary (Eg ) . Such a result corroborates the hypothesis made in Chatterjee et al. [17] Suggesting the introduction of an 

additional activation energy at the boundary. Furthermore, we notice that increasing the heating rate leads to decreasing 
activation energy (see Table 4). Indeed, by increasing temperature, the cages of zeolite channels undergo higher vibration and 
thermal agitation which allows the adsorbed hydrocarbon molecules to be diffuse more easily. 


Table 4 

Activation energy obtained from numerical program (kj.mol 1 ) 


Sample 

Rate (°C/s) 

Ethylene 

Toluene 

K 

E- 

T f (K) 

K 

E- 

T f (K) 

Na-ZSM-12 

3 

77 

129 

369 

84 

184 

416 

5 

70 

120 

378 

80 

180 

431 

9 

61 

110 

385 

75 

175 

438 

Ag-ZSM-12 

3 

(50*; 105”) 

(114*; 410**) 

391 

88 

187 

431 

5 

(43’; 73”) 

(111*; 388") 

410 

85 

185 

443 

9 

(34’; 55”) 

(107'; 350") 

413 

75 

180 

449 


Activation energy value obtained from the low T peak 
Activation energy value obtained from the high T peak 


For the Ag-ZSM-12 zeolite, numerical and experimental results for the desorption process are compared in Fig. 6. Since 
ethylene desorption curves (Fig. 6 (A)) exhibit two peaks corresponding, respectively, to Bronsted acid sites (low 
temperature peak) and Lewis acid centers (high temperature peak), a deconvolution operation was carried out using FitPeak 
Matlab code in order to uncouple peaks (Fig. 6). Fig. 5 (A and A’) show a benchmark between numerical and experimental 
results for the ethylene at 3°C/s as an example of a heating rate while Fig. 5 (B and B’) illustrate results for toluene. It is 
shown that numerical results are in a good agreement with experimental trends. The activation energy values used in 
Equations (9) and (11) for the Ag-ZSM-12 zeolite are reported in Table 4. The aforementioned remarks about the Na-ZSM- 
12 zeolite, namely: i) ethylene having activation energy lower than toluene, ii) activation energy inside the zeolite lower than 
the activation energy at the lattice boundary and iii) decreasing activation energy at increasing heating rate, still hold for the 
Ag-ZSM-12 zeolite. 
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Finally, it may be concluded that, for both ethylene and toluene, activation energies for Na-ZSM-12 are lower than those of 
Ag-ZSM-12 (see Table 4). This can be explained by the presence of strong Lewis acid sites in Ag + exchanged zeolite while 
the Na-ZSM-12 zeolite is characterized by the presence of mainly Bronsted acid sites. 







Fig.6. Ethylene thermodesorption profile of Ag-ZSM-12 with (A) 3°C/s, (B) 5°C/s and (C) 9°/s 

HEATING RATE 

V. Conclusion 

Diffusion of a single adsorbed component ethylene or toluene, inside the microporous ZSM-12 zeolite was studied 
experimentally and numerically. The Na + and Ag + forms of zeolite ZSM-12 were synthetized and characterized. The 
thermodesorption process was conducted at three different heating rates (3, 5 and 9°C/s) selected to be similar to the muffler 
heating rate during the cold start process. 

Experimentally, it was verified that the amount of both ethylene and toluene detected are essentially not affected by the 
heating rate. Contrary to the ethylene adsorbed on Ag-ZSM-12, where this molecule desorbed as two peaks, only one peak 
was observed for toluene. Moreover the amounts of ethylene adsorbed on ZSM-12 are low compared to those of toluene, due 
to the small size of this molecule in comparison to the pore diameter of this zeolite adsorbent. 

Numerically, using a transient diffusion boundary value problem coupled with Fick’s law, it was possible to explain all 
qualitative features of the temperature dependent desorption profiles observed experimentally. In this model, a coefficient C 

was included in the diffusion coefficient Co(t) (6) that has units of m 2 /s and consequently the term p[t , a) which describes 
the probability that a molecule occupies a site at distance x in the zeolite channels at time t is dimensionless. 

For the two solids considered in this study, activation energies E w and E e for toluene are much higher than those for 

ethylene owing to its bigger size. We have also shown that the values of these parameters are higher for the Ag + exchanged 
zeolite than for Na-ZSM-12 because of the presence of strong Lewis acid sites in the former material. In addition each 
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activation energy for Na-ZSM-12 and Ag-ZSM-12 decreases at increasing heating rate. This was explained by thermal 

vibrations of lattice which allow higher diffusivity at high temperature. 
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